Dynamical phase transitions after quenches in non-integrable models 



O 
(N 

D 

in 



C. Karrasch^'^ and D. Schuricht^ 
^Department of Physics, University of California, Berkeley, GA 95720, USA 
^Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA ana 
'^Institute for Theory of Statistical Physics, RWTH Aachen University and 
JARA-Fundamentals of Future Information Technology, 52056 Aachen, Germany 

We investigate the dynamics following sudden quenches across quantum critical points belonging 
to different universality classes. Specifically, we use matrix product state methods to study the 
quantum Ising chain in the presence of two additional terms which break integrability. We find that 
in all models the return probability to the initial state becomes a non-analytic function of time in 
the thermodynamic limit. This so-called 'dynamical phase transition' was first observed in a recent 
work by Heyl, Polkovnikov, and Kehrein [arXiv: 1 206. 2505 v2] for the exactly-solvable quantum Ising 
chain, which can be mapped to free fermions. Our results for 'interacting theories' indicate that 
non-analytic dynamics is a generic feature of sudden quenches across quantum critical points. We 
discuss potential connections to the dynamics of the order parameter. 
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I. INTRODUCTION 

One of the central quantities in statistical mechanicaii^ 
is the canonical partition function 



(1) 



where /? is the inverse temperature, L denotes the system 
size, and f{(3) is the free energy density. A phase transi- 
tion at a temperature /3c is defined as a non-analytic point 
of the free energy density. For real temperatures and fi- 
nite systems the partition function ([1]) is analytic, thus 
precluding the existence of a phase transition. In the 
thermodynamic limit, however, the free energy density 
may possess singularities. One way to characterize phase 
transitions was proposed by Yang and Lcci^ The starting 
pointii^ is the observation that the partition function will 
have zeros in the complex /3-plane. For a finite system 
these zeros are isolated and do not lie on the real axis. 
For L — OO, however, the zeros may coalesce into lines 
which can cut the real axis at the critical temperature 
Pc- This picture was established in detail for the two- 
dimensional Ising modelii^ as well as related systems. 

Recently, Heyl et al^ pointed out the formal similar- 
ity between the partition function ([T]) and the overlap 
amplitude 



G(t) = (*o|e- 



-iHt 



(2) 



For a quantum quench G{t) is the Loschmidt amplitudCfS. 
i.e. the overlap of the initial state 1^*0) with its time evo- 
lution e~'^* 1^0) under the post-quench Hamiltonian H . 
Specifically, Heyl et al. studied the analytic properties of 
the boundary partition function 



Z{z) = (*o|e 



(3) 



as a function of z <E C and established close analogies 
with equilibrium phase transitions in statistical mechan- 
ics that we outlined above. For real z = R, Eq. ([3]) can 



be interpreted^ii as the partition function of a system 
of length R with boundary conditions described by the 
boundary state |^'o); for z = it one recovers the overlap 
amplitude of Eq. ^ . 

Heyl et al^ specifically investigated the prototyp- 
ical quantum Ising chain which exhibits a quantum 
phase transition^ between ferromagnetic and paramag- 
netic ground states. The Ising chain can be mapped to 
a model of free fermions and thus allows an exact evalu- 
ation of the boundary partition functionj^i^ If |^'o) and 
H are associated with the same phase, the lines of zeros 
of Eq. (U) do not intersect the imaginary axis, while for 
a quench across the quantum phase transition zeros with 
z = it exist. The vanishing overlap of the time-evolved 
and the initial state translates into non-analyticities of 
the return amplitude Ll{t) = — ln|G'(i)p and can thus 
be viewed as a dynamical phase transition^ Similarly to 
its thermal counterpart, the non-analyticities in time only 
appear in the thermodynamic limit. They lead to a break- 
down of short-time expansions analogously to the break- 
down of high-temperature expansions in the vicinity of 
thermal critical points. The characteristic times at which 
the non-analyticities in the return amplitude occur were 
also observed in the dynamics of the order parameter 

In this paper we address the following question: Is a 
dynamical phase transition a generic feature of quantum 
quenches across quantum critical points? According to 
our numerical results, the answer is clearly yes. From 
the study of the quantum Ising chain this is not a priori 
clear, as the Ising chain can be mapped to an integrable 
model of free fermions whose dynamics might be special. 

To be more precise, we use the time-dependent density- 
matrix renormalization group (DMRG) to investigate 
quenches in the quantum Ising chain (where the com- 
parison with the exact solution allows to test our nu- 
merics), the transverse axial next- nearest- neighbour Ising 
(ANNNI) model (which corresponds to the Ising model 
complemented by additional next-nearest neighbor in- 
teractions), and a generalized quantum Ising chain in a 
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longitudinal magnetic field. While the quantum critical 
points in the first two models are of Ising type, in the 
latter model it belongs to a different universality class. 
Furthermore, the ANNNI model as well as the general- 
ized Ising chain cannot be mapped to a non-interacting 
theory and are, to the best of our knowledge, not inte- 
grable. In all three models we consistently observe that 
when quenching across the quantum critical point the re- 
turn amplitude shows non-analyticities at critical times 
t*, while for quenches within the same phase the time 
evolution of the return amplitude is completely smooth. 
Our results for 'interacting theories' thus indicate that 
non-analytic dynamics is a generic feature of quenches 
across critical points independently of their universality 
class. We also discuss potential connections to the dy- 
namics of the order parameter. 

This exposition is organized as follows: In Sec. |lT] we 
first review the results of Ref. for the quantum Ising 
chain. We test our DMRG numerics against the exact 
solution and present new data for a quench starting out 
of a spin-polarized ground state. In Sec. IIIII we discuss 
DMRG results for the ANNNI model. The behavior of 
the return amplitude is very similar to the quantum Ising 
chain - we observe non-analyticities even if the model is 
strongly non-integrable. In Sec. IIVI we study a general- 
ized quantum Ising chain whose quantum phase transi- 
tion does not belong to the Ising universality class. For 
this model we also find non-analytic behavior of the re- 
turn amplitude, although for certain parameters the non- 
analyticity may appear only at late times. We summarize 
our results in Sec. |Vl Technical details including the im- 
plementation of the DMRG calculations are elaborated 
on in the appendix. 

II. QUANTUM ISING CHAIN 

We begin our study with the quantum Ising chain 

Rising = -JY1 + ' 

i 

where we assume J > and g > 0- The model exhibits 
a quantum phase transition^ at = 1, which separates 
a ferromagnetic (FM) phase for g < gc from a paramag- 
netic (PM) phase for g > gc- In the FM the system pos- 
sesses two degenerate ground states |±) with (af) ^ 0, 
while the PM ground state with (ct|) = is unique. 
Close to the quantum critical point g = gc the correlation 
length ^ diverges as f ^ \g~ gc\~'^ with v = 1. Quenches 
across this quantum phase transition have been studied 
by various authorsj2r,iJ, In particular, Heyl et al^ inves- 
tigated the analytic properties of the boundary partition 
function ^ and showed that its lines of zeros cut the 
imaginary time axis z = \t for quenches across the quan- 
tum phase transition, giving rise to a dynamical phase 
transition at a critical time t* /2 [see Eq. below]. 

It is well-known that the quantum Ising chain can be 
mapped to a model of non-interacting fermions via a 



Jordan- Wigner transformation (see e.g. Ref. Hoh . and is 
thus certainly integrableJ^ In the fermionic languange, 
the Hamiltonian can be easily diagonalized, 

-ffising = ^eg(fc)ai.afc , (5) 

k 

where 



eg(fc) = 2jy (g- cosfc)2 -Hsin^fc , (6) 

and aj^ and are fermionic creation and annihilation 
operators. The fermions fulfill either anti-periodic or pe- 
riodic boundary conditions, usually referred to as Neveu- 
Schwarz (NS) or Ramond (R) sectors, respectively. The 
momenta k are correspondingly quantized as either half- 
integer or integer multiples of 27r/L. Clearly, at g = 
gc ~ 1 the dispersion in Eq. ([6]) becomes gapless, in- 
dicating the existence of the quantum phase transition. 
We note that the mapping from the spin model (U) to 
the fermionic model ([5]) involves a non-local transforma- 
tion of the spin operators, which leads to subtleties for 
quenches originating in the FM (see below). 

The quench protocol in the Ising model is implemented 
by suddenly switching the transverse field from its initial 
value g = go to its final value g = gi- The boundary par- 
tition function for the equivalent quench in the fermionic 
model can be calculated explicitly^i^ and, up to a con- 
stant, in the thermodynamic limit reads 

Z{z) = cxp[-L,nz)] , (7) 

with 

/W = - ln[cos2 0,+sin2 0fee-2-.iW] , (8) 

Jo ^TT 

where (pk = ^go (^) ~ ^gi (^) ^^'^ ^g{^) = arctan[sin fc/(g — 
cosfc)]/2 e [0,7r/2]. We stress that Eq. Q only applies 
to quenches starting from the unique ground state of the 
fermionic model, which lies in the NS sector for any finite 
system. In the PM phase of the Ising model this state 
corresponds to the unique ground state, while in the FM 
phase it corresponds to a superposition of the degener- 
ate ground states ji)*^ Hence, quenches originating in 
one of the spin-polarized states |±) are not described by 
Eq. (El). 

The zeros of Eq. ([7]) are located along the lines^ 
zn{k) = \ [ln(tan2 (jik) + i7r(2n + 1)] , n G Z , (9) 

where the argument of the logarithm in Eq. vanishes. 
For quenches across the quantum critical point these lines 
of zeros cross the imaginary time axis, that is there exist 
zeros with 5Hez„(fc*) = for k* = arccos[(l 4-170.91 )/(.9o + 
gi)]. This leads to a non-analyticity in the return ampli- 
tude 

l{t) = -y In \G{t)\^ = 2 mcf (it) (10) 
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infinite-system algorithm. Details of the numerical im- 
plementation can be found in App. |B] 

As discussed above, the ground state in the PM phase 
of the Ising model corresponds to the ground state of 
the fcrmionic model ([5]). Thus, Eq. ([8]) is applicable, 
and the return amplitude pO[) can be directly obtained. 
The comparison to the DMRG data is shown in Fig. [T] 
the agreement is good. In particular, we observe non- 
analyticities at the times t* defined in Eq. (fTTj) . In con- 
trast, for quenches within one phase the return amplitude 
is a smooth function without features at the times [we 
show an example for this behavior for the ANNNI model 
in Fig.Hb)]. 



FIG. 1. (Color online) Return amplitude l{t) defined in 
Eq. (fTU)) for a quench from the PM to the FM phase in the 
Ising model of Eq. Q. l{t) characterizes the overlap with the 
initial state in the thermodynamic limit. The DMRG data 
agree well with the analytic result obtained by mapping the 
model to free fermions [Eq. Q ; see also Ref. [3I . The return 
amplitude shows non-analytic behavior at the times indi- 
cated by the vertical dashed lines [Eq. . 

at the times 

t:-3mz„(r) = r + , (11) 

with 

t*- ^ - ^ / 90 + 91 ~ .^2) 

We note that either go < 1 and gi > 1 or vice 
versa. Eq. (fT2)) represents a new time scale generated 
by the quench. Eq. (|12p remains finite in the vicin- 
ity of the quantum critical point, where the system can 
be described^ in terms of massive Majorana fields (see 
App. 1^. This time scale was also observed^i^S as an os- 
cillation frequency in the time evolution of the order pa- 
rameter (erf) for quenches from the FM to the PM when 
starting in one of the polarized states |±). 

After this brief review we turn now to our DMRG cal- 
culations for quenches across the quantum critical point 
and show how the results discussed above can be recov- 
ered using this approach. 



A. DMRG results: Quench PM to FM 

In this section we discuss time-dependent DMRG data 
for quenches across the quantum critical point in the 
Ising model (|3]). We start with the simpler quench from 
the PM to the FM and discuss the more subtle quench 
from the FM to the PM in the next subsection. Unless 
stated otherwise, all DMRG calculations were carried out 
directly in the thermodynamic limit L — >■ 00 using an 



B. DMRG results: Quench FM to PM 

As a next step we consider a quench originating in the 
FM phase. We assume that the system is prepared in 
the state |+) with (erf) > and subsequently quench to 
the PM phase with gi > 1. The corresponding return 
amplitude is shown in Fig.[2l^a). We clearly observe non- 
analyticities. However, l{t) cannot be described using 
the analytic result of Eq. ([5]), and the non-analyticities 
no longer occur at the times t*. This can be attributed 
to the fact that the state | -I-) does not correspond to the 
ground state of the fermionic model (O; hence, Eq. ^ 
is not applicable. One can gain some analytic under- 
standing by considering a simple (but instructive) quench 
from 5o to (71 3> 1, where the time evolution after the 
quench is governed by H' = —Jgi J2i "'f (^cc also supple- 
mentary material to Ref. [^. The initial state for go = 
is simply given by |-|-) = [ft ■ • ■ t)i and a straightforward 
calculation yields 

;+(t) = -21n|cos(J9iO| . (13) 

Thus, non-analyticities occur at the times f„ = 2t*{n + 
1/2), t* = 7r/(2Jgi), which differ from the times i*„ for 
all m, n. Since the spins in H' do not interact at all, the 
return amplitude is perfectly periodic and does not decay. 
For general quenches starting from a spin-polarized state, 
such a periodicity cannot be observed [see Fig.[2Ja)]. 

The inset to Fig. [UJa) shows the time evolution of the 
order parameter (erf (t)). As reported previously^^i^ at 
sufficiently late times we observe oscillations with the 
frequency 2t*. For the simple quench discussed above 
we specifically obtain {<jf{t)) = cos{2Jgit). We note, 
however, that in the general case the relation between the 
times of non-analytic behavior in the return amplitude 
and the oscillation frequency of the order parameter is 
unclear (see the discussion for the ANNNI model below) . 

In order to make contact to the analytic result of 
Eq. also for quenches starting in the FM phase, we 
note that the ground state of the fermionic model ^ 
corresponds to a superposition of the FM ground states^ 
|NS) = {\+) - |-))/V2. DMRG data for a quench start- 
ing from this state is shown in Fig. [IJb) ; the agreement 
with the analytic result is again good. In particular, the 
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FIG. 2. (Color online) The same as in Fig. [T] but for quenches from the FM to the PM phase in the Ising model. In the 
thermodynamic limit, the ground state |±) within the FM phase is two- fold degenerate, (a) Quench performed starting from 
the polarized state |+). The return amplitude shows non-analytic behavior which, however, does not occur at the times 
defined in Eq. Inset: Time evolution of the order parameter {ai{t)), which at sufficiently late times oscillates with the 

frequency 21* . (b) Quench starting from the mixed state |NS) = (|-|-) — |— ))/\/2- The DMRG data agree well with the analytic 
result obtained for the corresponding quench in the fermionic model^ of Eq. ((5)). The return amplitude shows non-analytic 
behavior at the times . Inset: Comparison of the return amplitudes starting from the mixed and polarized states ]NS) and 
|+), respectively. 



return amplitude possesses non-analyticities at the times 
t* . As the initial state is an equal superposition of |+) 
and |— ), the order parameter (erf (i)) vanishes identically 
at all times. 

In the inset to Fig. E^b) we compare the return ampli- 
tude for quenches starting from the polarized state |+) 
and the mixed state |NS). We observe that for half of 
the time the return amplitudes are identical. This can 
again be understood by considering the simple quench 
introduced above. One finds that (if L is divisible by 
four) 



by the Hamiltonian 



Imit) 



■ ln[cos^{Jgit) + sin^(Jgit)] 



(14) 



For large systems this shows a switching behavior— de- 
pending on whether the first or the second term in the 
argument of the logarithm dominates, in complete anal- 
ogy with the DMRG data for generic quenches from the 
FM to the PM. The non-analyticities of Eq. ^ follow 
from cos^{Jgit) + sin^(Jgii) = 0, which in the thermo- 
dynamic limit L — J- oo yields t = . 

With this we conclude our analysis of quenches in the 
Ising model. In the remainder of this paper we address 
the question whether non-analytic behavior in the return 
amplitude can be observed for quenches across quantum 
critical points in other models. We begin by considering 
the ANNNI model in the following section. 



III. ANNNI MODEL 

As second model we investigate the transverse axial 
next-nearest-ncighbour Ising (ANNNI) modelii defined 



Ha 



NNNI 



= - J 



Aaf 



'i+2 



(15) 



Again we assume J > and g > 0, while A can be pos- 
itive or negative. Obviously, for A = we recover the 
quantum Ising chain of Eq. We note that Eq. ([15]) 
is invariant under g — > —g due to the transformation 
af'^ — > —af'^. Using a Jordan- Wigner transformation. 
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FIG. 3. (Color online) Sketch of the phase diagra m^^'^^ of 
the ANNNI model defined in Eq. (fTS)) as a function of A 
and g. There are four phases: a paramagnetic (PM) phase, a 
ferromagnetic (FM) phase, an anti phase (AP), and a floating 
phase (FP). The PM and FM phases are separated by an Ising 
transition located at (?c(A) defined in Eq. (|16p . We study 
quenches across the Ising transition as indicated by the solid 
arrow [see Fig.|4ja)], the dashed arrow [see Fig.|4jb)], and the 
dashed-dotted arrow (see Fig. [5]). 
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FIG. 4. (Color online) Return amplitude for a quench from the PM to the FM phase of the quantum Ising model in presence 
of integrability-breaking next-nearest neighbor interactions [the so-called ANNNI model; see Eq. (|15|l ]. (a) Quench in the 
transverse field g (indicated by the solid arrow in the phase diagram shown in Fig.[21). (b) Quench in the next-nearest neighbor 
interaction A (dashed arrow in Fig. [3]). In complete analogy with the integrable 'non- interacting' Ising chain (A = 0), the 
return amplitude exhibits non-analytic behavior as a function of time in the thermodynamic limit if one quenches across a 
critical point [note that the curve in (b) for Ai = 0.2 corresponds to a quench within the PM phase]. 



the ANNNI model can be mapped to a model of interact- 
ing fermions. To the best of our knowledge, the resulting 
system is not integrable and does not allow an exact so- 
lution like the quantum Ising chain. 

The phase diagram of the ANNNI model contains four 
phases (see Fig.[3])fi^— A paramagnetic phase (PM) with 
a unique ground state satisfying (erf) = 0; a ferromag- 
netic phase (FM) with doubly degenerate ground state 
with (erf) 7^ 0; an "anti phase" (AP) that schematically 
looks like tt||tt||; and a "floating phase" (FP) between 
the PM and the AP. The phase transition between the 
PM and the FM is in the Ising universality class with 
V = 1. For A < it is located at 



Fig. [5] shows DMRG data for the opposite quench from 
a polarized FM ground state to the PM phase (dashed- 
dotted arrow in Fig. . The appearance of kinks is again 
stable against interactions A 7^ 0. Even for the quan- 
tum Ising chain the kinks do not occur periodically if 
one starts from a spin-polarized state; this behavior be- 
comes more pronounced for A 7^ 0. In particular, the 
evolution between the kinks becomes highly non-trivial 
including smooth maxima and inflection points, suggest- 
ing that for such details interaction effects become im- 
portant and a simple picture based on the time evolution 
under a trivial Hamiltonian like H' is not sufficient to 
describe the dynamics. The order parameter, however, 



1 -h 2A = <7, 



2(H- A) 



(16) 



We will restrict ourselves to quenches across this phase 
transition in the following (see the arrows in Fig. [3]). 

We first concentrate on quenches from the PM to the 
FM phase. Fig. illustrates the effect of successively 
switching on the 'interaction' A for a quench analogous to 
the one shown in Fig. [1] (it corresponds to the solid arrow 
in Fig. [3]). The existence of non-analyticities is stable 
against interactions A. However, the critical times do 
not show any periodicity, i.e. it is not possible to write 
C = t*(A)(n -I- 1/2) with some interaction dependent 
time scale i*(A) replacing Eq. We have furthermore 
investigated an interaction quench with fixed go = gi = 
1.3 and Aq = 0, Ai > (depicted by the dashed arrow 
in Fig. [3]). The results are shown in Fig. UJ^b). For Ai = 
0.2, one does not leave the PM phase, and the return 
amplitude is a smooth function of time. In contrast, for 
Ai = 0.6 and Ai = 1 one enters the FM phase and l{t) 
becomes non-analytic as expected. Note that the model 
is strongly non-intcgrable for those parameters. 
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FIG. 5. (Color online) The same as in Fig.|4]but for a quench 
from the polarized ground state of the FM phase to the PM 
phase (dashed-dotted arrow in Fig. |3]). The return amplitude 
consistently features non-analyticities. Inset: The order pa- 
rameter oscillates periodically at sufficiently late times, but 
the precise connection of its dynamics to l(t) remains elusive. 
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FIG. 6. (Color online) Return amplitude for quenches within the generalized Ising model defined in Eq. (|17p . The model is 
critical at g — 1 for all (j) > 0. For > 0, the model has a unique ground state for any g, is not integrable, and the quantum 
critical point is no longer in the Ising universality class. One consistently observes non-analyticities in the time evolution; they 
occur periodically only at = 0. We note that for large angles <j) the non-analyticities appear at rather late times (see the 
inset) while the behavior for short times is absolutely smooth. 



still oscillates periodically at sufficieiitly late times (see 
the inset to Fig. [5]). A straightforward connection be- 
tween order-parameter dynamics and the appearance of 
non-analyticities in the return amplitude thus remains 
elusive. 

Our results for the ANNNI model show that the ap- 
pearance of a dynamical phase transition proposed in 
Ref.i is not an artefact of the quantum Ising chain: the 
non-analyticities in the time evolution of the return am- 
plitude are stable against integrability-breaking interac- 
tions. 



IV. ISING MODEL IN A MAGNETIC FIELD 

As last example we consider a generalized Ising chain 
in a longitudinal magnetic field i^^'^^ Specifically, we use 
the parametrization 

-ffgi = -J^[a,f(7.^i+< + (.9-l)(cos(/.<+sin(/.0] , 

(17) 

where J > 0, g > 0, and < \(j>\ < 7r/2. For </> = we 
recover the quantum Ising chain Eq. ([TT]) exhibits 
a quantum phase transition at g = 1. For any finite 
angle cj) the ground state is unique and, to the best of our 
knowledge, the model is not integrable. The transition 
at g = 1 is not in the Ising universality class; instead, the 
correlation length ^ diverges with the exponent i> — 8/15. 
The application of a Jordan- Wigncr transformation to 
the magnetic field term (g — 1) sine/) erf yields a non-local 
string of fermions. Thus, an analysis of Eq. (jl7p in terms 
of fermionic degrees of freedom is not possible. 

DMRG results for quenches at fixed (j>o = 4>i ~ (p from 
go > 1 to gi < 1 (which in the limit — )■ corresponds 
to quenches from the PM to the FM phase) and vice 
versa (polarized FM state to the PM phase) are shown 



in Fig.|6l The return amplitude features non-analyticities 
for arbitrarily large cj). As for the ANNNI model, these 
non-analyticities do not occur periodically and can po- 
tentially be shifted to rather late times. This is partic- 
ularly striking at = 7r/8 in Fig. IDJa) (note the inset). 
The dynamics of the generalized Ising model was also 
investigated by PoUmann et alr^ who studied a finite- 
velocity sweep across the quantum critical point instead 
of a sudden quench that we focus on. For the continuous 
sweep a non-intcgrable perturbation </) 7^ was found to 
lead to a broadening of the kinks. A better understand- 
ing of the difference between sudden and slow continuous 
quenches will be the topic of future investigations. 

Our results for the generalized Ising model, which is 
non-integrable and possesses a quantum phase transti- 
tion which is not in the Ising universality class, again 
indicate that the appearance of a dynamical phase tran- 
sition is a generic feature of sudden quenches across quan- 
tum critical points. 



V. CONCLUSIONS 

In this work we investigated sudden quenches across 
a quantum critical point of one-dimensional spin mod- 
els using the time-dependent density matrix renormal- 
ization group. We specifically studied the quantum Ising 
model in presence of two terms which break intcgrabil- 
ity (next-nearest neighbor interactions and a longitudi- 
nal field) and showed that in the thermodynamic limit 
the return probability Ll{t) = - In | (^ol e"'^* l^-o) P to 
the initial state |^'o) becomes a non-analytic function of 
time. This breakdown of short-time expansions is anal- 
ogous to the failure of high-temperature expansions of 
the partition function in the vicinity of equilibrium crit- 
ical points. It was thus termed 'dynamical phase tran- 
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sition' and investigated in detail for the integrablc (non- 
interacting) quantum Ising model in Ref. 0. The systems 
studied in our paper are not equivalent to free theories 
and not integrable. Thus, our results indicate that 'dy- 
namical phase transitions' are a generic feature of sudden 
quenches across a quantum critical points. For quenches 
originating in an ordered phase we studied the dynam- 
ics of the order parameter: it oscillates periodically at 
sufficiently late times, but a precise connection to the 
non-analyticities in the return amplitude remains elusive 
except for the quantum Ising chain. 
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Appendix A: Scaling limit in the Ising model 

For definitness we assume a quench from the PM to 
the FM phase, i.e. go > 1 and gi < 1. In order to study 
the scaling limit of Eq. (|12p we first introduce the masses 



Mo = 2J(.go - 1) , A/i = 2,7(1 - 51) , 



(Al) 



which are both positive. In terms of the masses the crit- 
ical time (|12p reads 



4 J + Mo - Ml 



iMo + Mi)Mi{4J - Ml) 



(A2) 



Now taking the scaling limit J — > 00, go,gi 1, and 
a — > (a is the lattice spacing previously set to one), 
while keeping the masses and the velocity v = 2Ja fixed, 
we find 



t* 



^/(Mo + Mi)Mi 



(A3) 



The existence of this time scale can also be inferred from 
the calculation of the boundary partition function Eq. ([3]) 
directly in the Ising field theory. For a quench across 
the quantum critical point the derivation requires the 
introduction^ of an ultra-violett cut-off in the form of 
an 'extrapolation time' as well as the regularization of 
infinite- volume divergencies fS^iS, 



We note that the time scale ()A3|) diverges as* t* oc 
1/y/Mi if one takes the conformal limit of the system 
after the quench Mi — > (i.e. one quenches to the quan- 
tum critical point), while quenches away from the quan- 
tum critical point (A/q ^ 0) simply lead to the finite 
value t* = tt/Mi. 
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FIG. 7. (Color online) DMRG calculation of the return am- 
plitude per site l{t) for a quench across the phase transition of 
the generalized Ising model of Eq. (|17|) with 00 = (/"i = 7r/32. 
As the system size L is increased one successively approaches 
the result obtained directly in the thermodynamic limit L = 
cxj via an infinite-system DMRG algorithmi^ 



Appendix B: DMRG calculations 

We study the quench dynamics of the quantum Ising 
model Eq. ^ and its non-integrable generalizations 
Eqs. ([TS]) and ^ using a DMRG algorithmJ^ The 
ground state of any one-dimensional system can be ex- 
pressed in terms of a matrix product state (MFS);^ 



\0102 ■ 



(Bl) 



where \a) denotes the local product basis characterized 
by . In the thermodynamic limit L — >■ 00 it is con- 
venient to directly work with a translationally-invariant 
statei^ 



(B2) 

The matrices ^[^■ ■^1'^^ are associated with a unit cell 
of size (e.g., A^ = 2 for the Ising model). We first 
determine the ground state l^^o) (corresponding to some 
initial parameters goj ^o, (/)o) by applying an imaginary 
time evolution exp(— ri/o) to a random state until the 
energy has converged to ten relevant digits (since we are 
only considering situations where Hq is gapped, this is 
very efficient; sec below). Thereafter, we carry out the 
real time evolution exp(— ii7it)|^'o) using a Hamiltonian 
Hi which features different parameters 51, Ai, 0i. 

In order to implement the above procedure, we em- 
ploy a standard time-dependent DMRG algorithm.— Af- 
ter factorizing the (real or imaginary) evolution operators 
exp(— cA-ff) using second or fourth order Trotter decom- 
position, they can be successively applied to Eqs. (jBl[) 
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or (jB2p . At each step AA, singular value decompositions 
are carried out to update (some of the) matrices A. We 
fix their dimension to a constant value x which consti- 
tutes our main numerical control paramater. The ground 
state of a gapped system only features finite entangle- 
ment and can thus be expressed exactly in terms of a 
MPS with finite x (we are always starting with a gapped 
Hq). During the real-time evolution using exp(— iiJit), 
the entanglement increases (we have carried out a global 
quench). We ensure that x is chosen large enough (i.e., 
the discarded weight is small enough) so that our results 
are numerically exact on the time scales that we study. 

For a finite system described by Eq. (jBl[) . the over- 
lap amplitude per site can be computed straightfor- 



wardly as the L-th root of the overlap between \^q) 
and exp(— iiJi^) |^'o). The overlap per site of two infi- 
nite MPS [described by matrices A and A in Eq. (jB2p ] is 
given by the iV-th root of the dominant eignvalue of the 
transfer matrix 

r({ai,ai},{ajv,aiv})= V (a™-^ . . . AI'^I-") * 

(7\...(7N 




(B3) 

A comparison between the results obtained by finite- and 
infinite-system algorithms is shown in Fig. [T] 
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